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Abstract 


In the past years, global navigation satellite systems (GNSS) have gained the core position 
concerning the geolocalization applications and services in urban environments. The 
major issue of the GNSS-based urban application expansion is related to the positioning 
service quality assurance, expressed in terms of accuracy, integrity, availability, and conti- 
nuity of the localization service. The dense urban environments, such as city centers, are 
challenging to the GNSS signal reception causing the frequent blockage of the line-of-sight 
(LOS) signals and the multipath phenomenon, referred to as the reception of the diff- 
racted/reflected echoes of the transmitted signal. These effects severely affect the pseudo- 
range and Doppler measurements, used by a GNSS receiver for the user’s position com- 
putation, which will further induce the computation of an erroneous positioning solution 
by the navigation processor down to a positioning loss in the presence of limited satellite 
visibility and few provided measurements. Therefore, advanced signal processing tech- 
niques do represent viable solutions aiming at the mitigation of these undesired effects in 
order to foster the accuracy and availability of the localization solution. This chapter will 
address in details the GNSS vector tracking (VT) receiver’s configuration able to cope with 
the urban environment-induced effects. 
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1. Introduction 


Global navigation satellite systems (GNSS) are increasingly present in our life and represent a 
key player in the world economy mostly due to the expansion of location-based services. An 
important part of the GNSS applications is found for the automotive usage in urban environ- 
ments that are characterized by difficult signal reception conditions. These services do exhibit by 
very stringent quality of service demands. In these harsh environments, the received signals are 
severely affected by the urban obstacles including buildings, road infrastructure, and foliage 
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generating strong signals’ fading and fast phase oscillations. Consequently, the generated 
pseudo-range and Doppler measurements are degraded. In the worst-case scenario, the direct 
signal can be totally blocked by the urban obstacles generating the GNSS signal blockage 
phenomena that introduce large biases on the pseudo-range measurements and cycle slips for 
the carrier phase observations. 


Standard GNSS receivers, which track each satellite independently through the technique 
referred to as scalar tracking, fail to cope with these harsh urban conditions. Therefore, this 
work is particularly focused on the proposal and detailed design of the dual-constellation 
single-frequency vector tracking architecture for the automotive usage in urban environment. 
This architecture can improve the tracking of some attenuated or blocked signals due to the 
channel-aiding property based on the navigation solution estimation and thus assure a robust 
navigation solution estimation. 


2. GNSS receiver structure 


This section aims at providing a detailed description of the GNSS receiver structure with an 
emphasis on the GNSS signal processing stage from the GNSS signal reception up to the 
measurement generation process. 


2.1. GNSS digital signal processing 


The digital signal processing is the consecutive stage of the analog front end, taking as input of 
the sampled and discretized received signals from each satellite in view. Herein, the digitized 
signal is fed to multiple channel processing blocks, corresponding to each received satellite 
signal. The digital signal processing blocks conduct three main operations, the correlation, the 
acquisition, and the code/carrier tracking [1]. The scalar tracking process is achieved in a 
channelized structure for all the satellites in view. The objective of the tracking process is to 
refine the coarse estimations of the code delay and Doppler frequency provided by the acqui- 
sition block and to precisely follow the signal properties changing over time [2]. The GNSS 
scalar tracking architecture is illustrated in Figure 1. 
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Digitized P, Op 


signal 


Figure 1. The conventional GNSS tracking architecture. 
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The scalar tracking process, employed for each tracked satellite, comprised the following main 
modules: 


e Code tracking: in charge of estimating the code delay offset (€+) between the incoming 
signal’s code and the locally generated replica based on the closed feedback loop referred 
to as the delay lock loop (DLL). In this process, at least three delayed code replicas, known 
as the early, prompt, and late replica, are generated according to the correlation process 
with the incoming signal. 


e Carrier tracking: responsible of estimating the residual Doppler shift (ep) and the carrier 
phase offset (£p). The carrier tracking module that provides the estimation of the Doppler 
shift (ep) is called as frequency lock loop (FLL), while the carrier phase error (£p) com- 
pensation is conducted by the phase lock loops (PLL) [2]. It is important to highlight the 
fact that for certain applications with important user’s dynamics, an FLL-aided PLL may 
be also employed. 


The scalar tracking architecture, whose high-level representation is provided in Figure 1, 
includes the following main processing blocks: 


e = Correlators: accumulating the combination of the sampled signal with the code and carrier 
replicas through mixing the in-phase and quadrature signal branches with the three delayed 
(early, prompt, and late) code-spreading sequences from the GNSS code generator block 


e Codelcarrier discriminators: in charge of estimating the code delay and carrier frequency/ 
phase errors based on the three correlator output pairs 


e Low-pass filters: filtering the discriminators’ outputs for noise reduction at the entry of the 
code/carrier local oscillators 


e = Numerical control oscillator (NCO): operating in a feedback loop manner and providing the 
required correction to the code and carrier generators based on the filtered discriminator 
output, which is further used to generate the local replicas for the successive measurement 
epoch 


The association of the low-pass filters and the NCOs provides the equivalent loop filter, whose 
response to the user’s dynamics is strictly related to two parameters such as the filter’s loop 
order and the equivalent noise bandwidth. In details, a higher noise bandwidth implies faster 
loop response time (and thus better response to high user’s dynamics) but with the drawback 
of dealing with noisier results due to the shorter integration time. Considering the filter’s loop 
order, the higher the loop filter, the better the filter capability to follow the high-order user’s 
dynamics is [3]. Keen readers may find detailed description of the loop filters in [3-5]. 


The scalar tracking is continuously run for each satellite-user channel in order to precisely 
estimate the code delay and carrier frequency/phase evolution in time [2]. In the GNSS 
receiver, the code delay and carrier frequency/phase lock loops are jointly used. 


2.1.1. GNSS code delay tracking 


The code delay tracking process is directly performed after the acquisition stage, responsible 
among other tasks of the detection of the incoming signal. The code delay tracking process 
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assures the alignment between the incoming signal’s code and the local replica’s code by 
refining the code delay error measurement from the discriminator block that further steers the 
code NCO [5]. The code tracking is performed by means of a DLL that is a feedback loop 
capable of steering the local PRN code delay based on the estimation of the code delay error £, 
[6]. This code delay error measurement is later used to compute the pseudo-range observation. 


The correlation between each local replica with the in-phase and quadrature signal samples 
generates one correlator pair. Finally, three correlator pairs are obtained at the end of this 
operation, expressed by 


A keTe 
IE, = z Ae Re (ex k+ 5 ). cos (ëp; k) sinc(7 Efnt) + NIE,k 


A ; 
IP, = z Ae Rele) cos (Ek) . sinc( rep, t) + NIP, k 
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where 


e The triplet (ex; Eg ks Efp D) represents the code delay, carrier phase, and frequency esti- 
mation errors at the current epoch k. 


e k.-T, refers to the early-late code chip spacing with k, and T, representing the fraction of 
chip spacing and the code chip period, respectively. 


e — ny denotes the noise term at the correlator output (where x stands for the in-phase (I) or 
quadrature (Q) and y for the early (E), prompt (P), or late (L) code delays), which is 
correlated and Gaussian distributed with zero mean and variance i 

The most common discriminators employed in GNSS receivers are the noncoherent early- 

minus-late power (EMLP) and the dot product (DP) discriminators due to their tracking 

robustness since insensitive toward the carrier phase information. These two discriminators 


are defined by 


Demup(€x,k) = (IE + QE?) — (IL? + QLK) 


2 
Dpp(€r,k) = (IEx — ILx)-IPx + (QEx — QLx)-QP; a 


The chip spacing is a crucial parameter for the discriminator function and thus must be 
selected so that the early and late correlator outputs are always evaluated at the correlation 
function main peak. Therefore, the EMLP and DP discriminators require a correlator spacing 
set less than 1 chip and 0.5 chip for the GPS L1 binary phase-shift keying (BPSK) and Galileo 
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binary offset carrier (BOC) signals, respectively [3]. In order to remove the amplitude sensitiv- 
ity of the code discriminators, normalization factors shall be applied to have a direct access to 
the code tracking error according to [7, 8]. 


2.1.2. GNSS carrier phase tracking 


The carrier phase tracking is achieved by the phase lock loop (PLL) that is responsible of 
keeping the carrier phase alignment between the incoming signal and the locally generated 
replica. The main goals of the phase tracking loop are: 


e The computation of a phase reference for the detection of the GNSS modulated data signal 
e The provision of precise Doppler measurements by using the phase rate information 


To fulfill these objectives, the PLL uses a carrier phase discriminator to estimate the resultant phase 
estimation error (€,(k)) between the incoming signal phase and the replica phase, according to the 
in-phase and quadrature prompt correlator output pairs (IP, QP). Afterward, the PLL filter filters 
out the noise, and afterward, the carrier NCO transforms the estimated phase error into a 
frequency offset that modifies the NCO nominal frequency for the successive epoch [6]. 


The main phase discriminators used for the data channels are the following [4]: 


e Dot product (DP) or the generic Costas discriminator 
Dpp (£p,k) = QP IP (3) 
e The extended arctangent (Atan2) phase discriminator [9] 


P 
DAtan2(€p,k) = atan2 ($) (4) 


3. Vector tracking architecture 


3.1. Urban environment-induced effects 


The urban environment causes important threats to the GNSS signal reception and its posteriori 
processing, severely deteriorating the navigation solution accuracy and availability. The prob- 
lems of the urban environment in the navigation domain can be summarized as follows: 


e Multipath: denoted as the joint reception of both the direct LOS signal and the reflected or 
diffracted GNSS LOS echoes from the urban obstacles 


e Attenuation or blockage of the GNSS LOS signal: arising due to the partial or total obstruction 
of the GNSS LOS from the urban environment characteristics 
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e Interference: resulting from the presence of wide signal sources transmitting in the adjacent 
GNSS frequency bands 


The resulting consequences of these urban environment error sources from the GNSS receiver’s 
point of view are the following [10]: 


e Distortion of the correlation function: computed between the received multipath-contaminated 
signal and the NCO-generated replica. The tracking of the multipath-affected signals leads 
to a deterioration of code delay and carrier frequency/phase estimation accuracy down to 
the loss of lock of the tracking loops. As a result, the estimated pseudo-ranges and pseudo- 
range rate measurements from the code and carrier tracking loops, respectively, which are 
fed to the navigation filter are significantly corrupted. 


e = Only non-LOS (NLOS) signal reception: resulting from the total obstruction of the direct 
LOS GNSS signals and thus inducing large biases on the pseudo-range and Doppler 
measurements if only NLOS signals are tracked. 


Finally, the resulting degraded measurements cause the navigation processor to calculate an 
inaccurate position solution or even to be unable to compute one in the case of few available 
measurements. 


3.2. Vector tracking motivation 


Vector tracking algorithms represent advanced GNSS signal processing methods, having the 
ability to function at lower carrier-to-noise power (C/No) ratios and in higher user’s dynamics 
than traditional GNSS receivers [11]. Contrary to the conventional or scalar tracking architecture, 
where each visible satellite channel is being independently tracked, vector tracking performs a 
joint signal tracking of all the available satellites. It exploits the knowledge of the estimated 
receiver’s position and velocity to control the tracking loops’ feedback [6]. The comparison 
between the scalar tracking and vector tracking architectures is illustrated in Figure 2. 


Concerning the vector tracking architecture, depicted in Figure 2(b), the code and carrier 
tracking loops of all the visible satellites are connected via the navigation solution computed 
by the navigation filter. The individual code/carrier loop filters and NCOs, illustrated by the 
dashed red line in the left figure, are abolished and substituted by the vectorized code/carrier 
update block depicted in blue. 


3.3. Proposed vectorized GNSS receiver architecture 


Aiming at a robust GNSS signal tracking and navigation technique in urban environment, a 
dual-constellation GPS + Galileo single-frequency L1/E1 vector delay/frequency-locked loop 
(VDFLL) architecture is proposed and implemented for the automotive usage, whose architec- 
ture is illustrated in Figure 3. 


In fact, the use of dual-constellation measurements significantly improves the availability of a 
navigation solution in urban canyons. Moreover, the implementation of the VDFLL tracking 
architecture, where the navigation filter is responsible of estimating both the code delay and 
the Doppler frequency change of each received signal, enhances the vehicle dynamic tracking 


Robust GNSS Positioning in Urban Environment 
http://dx.doi.org/10.5772/intechopen.80412 


capability of the GNSS receiver [6]. As it can be clearly observed in Figure 3, the central 


navigation filter accepts the code (e0) and carrier (4?) discriminator outputs for each GPS 


(i = 1 + N1) and Galileo (i = 1 + N2) tracked channel as its input vector. Furthermore, the code 
and carrier NCO update process is achieved by projecting the predicted navigation solution 
into the pseudo-range and pseudo-range rate domains. 


The VDFLL navigation block employs a Kalman filter (KF) that is a Bayesian estimation 
technique, which incorporates the measurements from the past epochs to obtain a more 
accurate navigation solution. In this work, the selection criteria of the extended Kalman filter 
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Figure 2. The high-level comparison of (a) conventional or scalar tracking and (b) vector tracking architectures. 


Figure 3. The noncoherent L1/E1 VDFLL architecture. 


35 


36 Accuracy of GNSS Methods 


(EKF) algorithm are related to its capability to resolve the nonlinearity issues for the GNSS 


navigation system. 


The EKF estimation process is composed of two main stages: 


e State prediction: performing the time propagation between two consecutive epochs (k — 1) 


to the current one k of the state vector X,,_; and its covariance matrix Pxx-1 


e Measurement correction: refining the a priori state vector and covariance matrix estimations 
(Xige-1,Pke-1) by feeding the current epoch measurements (Zjnput) into the filter and thus 


obtaining the improved a posteriori estimates (Xiks Prk) 


3.3.1. VDFLL EKF state space description 


The continuous-time EKF state model is given by 


d 
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where 
e X is the state vector with the following entries to be estimated such as the user’s position 
[x y z], velocity [xyz], and user’s clock bias and drift terms [brb Re] expressed in unit of [m] 
and I, respectively. 
e F denotes the state transition matrix characterizing the user’s motion and clock dynamics. 
e B represents the colored noise transition matrix. 


e w is the process noise vector reflecting the user’s dynamics (position and velocity) and 
receiver’s oscillator (clock bias and drift) uncertainties affecting the system model 


0 
A2x2 = | 


0 A represents the position/velocity and clock biases/drift state transition sub- 


matrices. 


The main tuning factors of the process noise vector w and its associated covariance matrix Q 
can be summarized into two main categories according to their nature: 


e User's dynamics: expressing the vehicle's dynamic uncertainties concerning the velocity 


error variance terms along the three ECEF axes (<2 , oy o | that are further projected in 
the position domain through the state transition sub-matrix A2x2 


e  Receiver’s oscillator noise: expressed in terms of the oscillator’s phase and frequency noise 
spectral densities affecting the receiver’s clock o? and drift 07 biases 


Passing to the discrete time domain, the system or dynamic model of the VDFLL navigation 
filter can be detailed as follows: 


Xp = D- Xka + Wy (6) 


where X; is the state vector forward projection from the k — 1" to the k" time epoch and ® 
represents the dynamics of the user platform and clock, expressed as follows 


@=I+F-AT, (7) 


where AT = ti_1 — t is the time step between two successive epochs. 


Substituting Eq. (7) into the continuous state transition matrix of Eq. (5), the final discrete state 
transition matrix is provided by 


Adg2x2 2x2 02x2 02x2 
D= 02x2 Aa2x2 02x2 02x2 (8) 
022 O2x2 Age2x2 O2x2 


02x2  Oox2 02x2 Agoxa 
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1 H 


where Aq 2x2 = p 1 


The process noise covariance matrix Q, = diag [Q b Qy, k Qz k» Qerk, A in the discrete domain 


per each entry can be expressed as 


Q; = E [wrw] 
=f Plte, T QC) DT (ty, t)-dt 


~ Jika 


(9) 


Therefore, the process noise discretization for the position and velocity states along the three 
reference axes is computed as 


tk 


Quiyizk = | Ad,2x2(tk, T)Qox2(T) Ad lt T)-dt 


tea 


f 1 e Uo -M | 1 i 
= : dt 10 
kwl 4 0 ome | LAT 1 (y 
AT? AT’ 
_ 2... /3 l2 
xlylz | ar ie AT 
where (o l a ; o; ) are the velocity noise error variances along the three navigation axes. 


Applying the discretization process of Eq. (9) to the user’s clock covariance states, the follow- 
ing relation is obtained: 


tk 


Quik k = | Aa2x2(tk, T)-Qatk, 2x2(T)Ag alik T)-dt 
tk—1 
i 1 AT] To i 1 ' 
<A . . ‘dT 11 
fo 14 e] Lan a (1) 
Pe, fs o7-AT? ") 


7 2AT? 2, 
og^ /2 og AT 


The discrete receiver’s clock process noise covariance matrix is modeled based on the Allan 
variance parameters [12]. 


3.3.2. VDFLL EKF observation model 


The nonlinear relation between the state and measurement vector in an EKF is provided by 


Zk = h(Xx) + Uk (12) 
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where 


e his the nonlinear function relating the measurement zz to the state Xx. 


e =v is the measurement noise vector that is modeled as a zero-mean uncorrelated Gaussian 
noise process and independent to the process noise wr. 


The measurement vector zy comprises pseudo-ranges p® and Doppler measurements 6") and 
output from the code/carrier tracking process for the i= 1 + N GPS L1/Galileo F1 tracking 
channels [12]: 


~~ (ep -p™) : (Dpp) (13) 


The GNSS pseudo-range measurements of a given satellite i (from the GPS (N1) and Galileo 
(N2) satellites in view) at epoch k are rewritten as 


, Iri(k) — ru(k)| + Xk(7) + e®cps(k) 0<i<N1 
p(k) = | (14) 
bri(k) — ru(K)| + Xe(7) + Oc alk), = NI <isN 


where 


; 2 ; 2 . 2 

© irk) —r4(k)| = i (a (k) — X(1)) $ GP (k) — X:(3)) $ (2? (k) — X,(5)) is the actual 
satellite-to-user Euclidian distance at the current epoch k where the triplets 
(20°, 20°) (k) and (X;(1), X.(3), X;(5)) represent the i" satellite and user’s coordinates 
in the ECEF reference frame, respectively. 

e X;(7) represents the receiver’s clock bias term w.r.t the GPS time. 

° ee sca lK) represents the error due to the receiver’s thermal noise, assumed to be white, 
centered Gaussian distributed. 


Meanwhile, the remaining N-entries of the measurement vector z;, constituted by the Doppler 
measurements from both the GPS and Galileo satellites, are related to the state vector through 


the observation function a) [6]: 


ny) (Xe) = (PE = X2) ) al? E) + (GP) = Xea) af 
| (15) 
+ (2, (&) = Xx(6)) al? + Xe(8) + E epsyca lK) 


where 
o (a\? ' ay? ; ak) (k) denotes the LOS projections along the three navigation axes. 


e  X,(8) is the receiver’s clock drift state common for both the GPS and Galileo Doppler 
measurements. 
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° a /Gal (k) denotes the receiver’s thermal noise effect on the GPS and Galileo pseudo- 


range rate measurements. 


The measurement noise vector v; is modeled as a zero-mean uncorrelated Gaussian noise 
process with the measurement noise covariance matrix Rx that has the following entries in the 
main diagonal: 


A . 
° Code_Discri, j for j=1-N 
Ry = wa (16) 


2 ; = ee 
© Carr—Discri, j for J= 1---N 


where the first entry refers to the pseudo-range error variance terms for the tracked GPS and 
Galileo satellites, while the second one is a common term for the pseudo-range rate error 
variance for all tracked satellites [12]. 


4. Description of the test setup 


Within the scope of this research, a realistic dual-constellation dual-frequency GNSS signal 
emulator, capable of simulating the GNSS signal reception at the correlator output level and 
including the navigation module, has been developed. The developed signal emulator is a 
powerful tool for flexible and reliable GNSS receiver testing, for which all the processing blocks 
from the GNSS signals’ correlation function, passing through the channels’ tracking module 
and up to the different navigation algorithms, are all designed in a modular manner [12]. 


Two different GNSS receiver architectures are herein analyzed with the scope of performance 
comparison: 


e The proposed L1/E1 VDFLL EKF architecture working at AT = 20 ms integration time and 
thus providing 50 Hz code and carrier frequency updates 


e The classic scalar tracking architecture employing a third-order loop PLL and a DLL, with 
an EKF positioning module for the navigation solution computation operating at the same 
rate as for the VDFLL EKF case (50 Hz) 


4.1. Test scenario 


The urban navigation test is herein performed on a car trajectory in Toulouse city center based 
on the collected data from the real test campaign in Toulouse urban area, by employing a 
NovAtel’s ProPak receiver mounted on the car. The recorded trajectory of 600 s duration is 
presented in Figure 4. 


The simulated signal reception conditions are that of a complete urban multipath model 
included to the receiver processing blocks with a maximum of 13 simultaneously tracked GPS 
L1 and Galileo E1 channels during the whole trajectory. 


The receiver tracking parameters employed in the test scenario, defining the scalar (ST) and the 
vector tracking (VT) algorithms, are summarized in Table 1. 
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Figure 4. The reference car trajectory. 


The navigation performance assessment of the proposed vectorized architecture with respect 
to the classic receiver configuration is performed in degraded signal reception conditions. 


4.2. Urban multipath model 


In this work, a widely used urban propagation channel model has been used to generate a 
representative of urban environment signal’s reception conditions. This model, known as the 
DLR land mobile multipath channel model, was developed thanks to an extensive measure- 
ment campaign conducted by DLR in Munich urban and suburban areas in 2002. This model is 


Parameters ST VT 


L1/E1 code delay tracking 


DLL order Ist Not used 

DLL configuration Carrier-aided DLL Not used 

GPS L1 chip spacing (kc_11) (chip) 0.5 

GAL E1 chip spacing (kc_r1) (chip) 0.2 

Discriminator Early-minus-late power (EMLP) 

DLL update period (s) 0.02 Position update rate 
DLL noise bandwidth (BprL-n) (Hz) 1 Not used 


Carrier phase/frequency tracking 


Carrier estimation Phase Frequency 

PLL order 3 Not used 

PLL update period (s) 0.02 Velocity update rate 
PLL noise bandwidth (BprL-n) (Hz) 10 Not used 
Discriminator type Atan2 Cross product (CP) 


Table 1. Tracking loop parameters. 
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a wideband propagation channel model where each LOS and multipath echo are individually 
considered [13]. 


Moreover, the DLR model is a hybrid statistic/deterministic mathematical propagation channel 
model. The statistical part refers to the generation of a random urban scenario from a given set 
of channel-defining parameters, as depicted in Figure 5. Afterward, the LOS and NLOS echo 
information, provided as the DLR model output, is fed in the tracking stage at the correlator 
output level per tracked satellite. 


4.3. Performance comparison results 


The position error comparison between the scalar tracking receiver and the VDFLL algorithm 
is presented in Figure 6, presenting the EKF estimation errors along the entire trajectory in the 
navigation frame expressed in the along- and cross-track coordinates, in blue the VDFLL and 
in red the KF with scalar tracking. Moreover, the blue- and red-dotted curves denote the (95%) 
or two o covariance bounds, where o is the estimation error standard deviation estimated by 
the Kalman filters for the two architectures, respectively. 


The position error plots in the vehicle frame (along- and cross-track coordinates), illustrated in 
Figure 6(a) and (b), demonstrate a clear superiority of the VDFLL algorithm, observed in terms 
of the low-position error variations along the trajectory that is an evident indicator of the 
VDELL capability in coping with the multipath conditions. Larger position errors are observed 
for both architectures in the cross-track coordinate in Figure 6(b) that is related to the higher 
multipath-induced bias on the lateral projection. A significant position bias of approximately 
30 m is seen for the scalar tracking receiver (in red) in the cross-track coordinate from the 70th 
to the 100th epoch, which coincides with strong satellite outage event and thus limited satellite 


~ \ Neen NNN 


Figure 5. Artificial urban scenario generated by the DLR urban propagation channel model [13]. 
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Along track position estimation error comparison 
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Cross track position estimation error comparison 
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Figure 6. Position error comparison between the scalar tracking receiver (in red) and the VDFLL algorithm (in blue) for 
the (a) along-track error [m] and (b) cross-track error [m]. 


measurements. Furthermore, the ST + EKF covariance bounds are significantly increased 
during this period due to the higher position estimation uncertainty since only four “good” 
measurements from the four locked LOS satellites are used for the navigation solution compu- 
tation. On the contrary, the proposed VDFLL algorithm assures a positioning stability and 
much tighter confidence bounds related to the interchannel-aiding capability exhibited by the 
VDFLL EKF filter. Finally, it can be seen that the vector architecture better assesses the confi- 
dence on the computed position. 
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5. Conclusions 


The chapter was written in the context of the GNSS use in urban environment that is particu- 
larly challenging to the GNSS signal reception due to multipath and direct signal blockages, 
which significantly affect the signal processing and further degrade the position accuracy and 
availability. For this matter, this work was focused on the design of a dual-constellation GPS/ 
Galileo and single-frequency L1/E1 band vector tracking architecture for automotive usage in 
urban environment. This architecture ensures a better receiver’s dynamic estimation due to the 
joint code delay and Doppler carrier frequency tracking for all the satellites in view performed 
by the common navigation EKF filter. The detailed flowchart of the proposed vector tracking 
algorithm and the relation between the state vector and observation model were also exposed. 


The generation of the urban environment multipath conditions was detailed in the second part of 
Section 4. In this work, the wideband propagation model referred to as the DLR land mobile 
multipath channel model was employed. Herein, the urban environment conditions were gener- 
ated separately for each GPS- and Galileo-tracked satellite by feeding their elevation/azimuth 
angles and the reference car trajectory to the DLR urban channel. The generated channel model 
samples were stored and directly fed to the signal emulator at the correlator level. 


Furthermore, the dynamic car trajectory along with the test parameters related to the scalar 
(ST) and vector tracking (VT) loop design was introduced. The emphasis of this work was 
dedicated to the performance analysis of the two architectures under study in the navigation 
level, expressed by the position estimation errors in the vehicle navigation frame. The VDFLL 
superiority was observed in the navigation domain, especially when referring to the position 
estimation accuracies. Nearly twice lower position estimation error variations were observed 
for the VDFLL architecture w.r.t the scalar receiver configuration. Furthermore, the VDFLL 
robustness and reactivity were noted especially during the satellite outage intervals that 
induce a reduced number of measurements fed to the navigation processor. During these 
intervals, stable and accurate navigation solution estimations are assured by the VDFLL filter 
thanks to the code/carrier NCO updates steered from the navigation solution estimations. 
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Acronyms and abbreviations 


ADC analog-to-digital converter 


AGC automated gate control 
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DLL delay lock loop 

DP dot product 

E early 

EMLP early-minus-late power 
EKF extended Kalman filter 
FLL frequency lock loop 

GNSS global navigation satellite systems 
I in phase 

IF intermediate frequency 

KF Kalman filter 

L late 

LNA low-noise amplifier 

LS least squares 

NCO numerical control oscillator 
P prompt 

PLL phase lock loop 

PVT position velocity and time 
Q quadrature phase 

SiS signal in space 

ST scalar tracking 

VT vector tracking 

VDFLL vector delay/frequency-locked loop 
WLS weighted least squares 
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